Magnetic correlations and quantum criticality in the insulating antiferromagnetic, 
insulating spin liquid, renormalized Fermi liquid, and metallic antiferromagnetic 

phases of the Mott system V 2 3 

Wei BacQ 

Physics Department, Brookhaven National Laboratory, Upton, New York 11973 

C. Broholm 

Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland 21218 
and Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 

G. Aeppli 

NEC, 4 Independence Way, Princeton, New Jersey 08540 
S. A. Carter 

Department of Physics, University of California, Santa Cruz, California 95064 

P. Dai 

Solid State Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 31831 

T. F. Rosenbaum 

James Franck Institute and Department of Physics, University of Chicago, Chicago, Illinois 60637 

J. M. Honig and P. Metcalf 
Department of Chemistry, Purdue University, West Lafayette, Indiana 47907 

S. F. Trevino 

United States Army Research Laboratory, Adelphi, Maryland 20783 
and Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 

(February 1, 2008) 

Magnetic correlations in all four phases of pure and doped vanadium sesquioxide (V2O3) have 
been examined by magnetic thermal neutron scattering. Specifically, we have studied the anti- 
ferromagnetic and paramagnetic phases of metallic V2- y 03, the antiferromagnetic insulating and 
paramagnetic metallic phases of stoichiometric V2O3, and the antiferromagnetic and paramagnetic 
phases of insulating Vi.944Cro.o5603. While the antiferromagnetic insulator can be accounted for 
by a localized Heisenberg spin model, the long range order in the antiferromagnetic metal is an in- 
commensurate spin-density-wave, resulting from a Fermi surface nesting instability. Spin dynamics 
in the strongly correlated metal are dominated by spin fluctuations with a "single lobe" spectrum 
in the Stoner electron- hole continuum. Furthermore, our results in metallic V2O3 represent an un- 
precedentedly complete characterization of the spin fluctuations near a metallic quantum critical 
point, and provide quantitative support for the self-consistent renormalization theory for itinerant 
antiferromagnets in the small moment limit. Dynamic magnetic correlations for huj < ksT in the 
paramagnetic insulator carry substantial magnetic spectral weight. However, they are extremely 
short-ranged, extending only to the nearest neighbors. The phase transition to the antiferromag- 
netic insulator, from the paramagnetic metal and the paramagnetic insulator, introduces a sudden 
switching of magnetic correlations to a different spatial periodicity which indicates a sudden change 
in the underlying spin Hamiltonian. To describe this phase transition and also the unusual short 
range order in the paramagnetic state, it seems necessary to take into account the orbital degrees 
of freedom associated with the degenerate d-orbitals at the Fermi level in V2O3. 

PACS number(s): 71.27.+a, 75.40. Gb, 75.30.Fv, 71.30. +h 



I. INTRODUCTION 

The Mott metal-insulator transition |l| is a localiza- 
tion phenomenon driven by Coulomb repulsion between 



electrons. Mott insulators are common among transition 
metal oxides and are an important group of materials 
which cannot be accounted for by conventional band the- 
ories of solids. Using a simplified model which now bears 
his name, 0] Hubbard demonstrated a metal-insulator 
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transition, in which a half-filled band is split by elec- 
tronic correlations into upper (unfilled) and lower (filled) 
Hubbard bands. 0] His approach, however, failed to pro- 
duce a Fermi surface on the metallic side of the Mott 
transition jf|. Also, there were unresolved questions con- 
cerning electronic spectral weight of the lower Hubbard 
band. Q Starting from the metallic side, a Fermi liq- 
uid description of the correlated metal was provided by 
Brinkman and Rice || . They found that while the Pauli 
spin susceptibility, x, and the Sommerfeld constant, 7, 
are both strongly enhanced by Coulomb repulsion, the 
Wilson ratio, xlli remains constant. Such behavior re- 
flects a strongly enhanced effective mass and it was found 
that the metal to insulator transition could be associated 
with a divergence in the effective mass of these Fermionic 
quasi-particles. 

The discovery of high temperature superconductivity 
in cuprates has revived interest in Mott systems Q. To 
be precise, the cuprates are classified as charge transfer 
systems [||, but similar low energy physics is expected 
for both classes of materials. Progress in experimental 
techniques, such as photo-emission and X-ray absorp- 
tion spectroscopy, has recently enabled direct measure- 
ments of the Hubbard bands and the Brinkman- Rice res- 
onance. U O n the theoretical side, there has also been 
progress especially in using infinite dimensional mean 
field theories [^0) to calculate the local spectral function 
of the Hubbard model close to the metal-insulator tran- 
sition [ fi"l| . A particularly important result is the synthe- 
sis of the seemingly contradictory pictures of Hubbard, 
Brinkman-Rice, and Slater on the metal-insulator tran- 
sition in strongly correlated electron systems fl^| . It is 
encouraging that these D=oo theories to a large extent 
reproduce the experimentally observed features [^| in the 
electronic density of states for three dimensional (3-D) 
transition metal oxides. 

While there has been a remarkable convergence of ex- 
periments and theories which probe the electronic density 
of states across the Mott metal-insulator transition, lit- 
tle work has been done to explore magnetism close to 
this phase boundary for materials apart from the lam- 
inar cuprates. On the insulating side, the Heisenberg 
spin Hamiltonian was expected to offer a good descrip- 
tion of magnetism but as we shall see in the following, 
this turns out not to be the case for systems with degen- 
erate atomic orbitals. On the metallic side of the tran- 
sition, one might have hoped that the local spin picture 
would survive because charge carriers in the Brinkman- 
Rice liquid are heavy and nearly localized. As we shall 
see this is also inconsistent with experiments. It turns 
out that the mobile quasi-particles have profound effects 
on magnetism in the metal. Indeed, the magnetism is 
that which one would associate with a quantum critical 
state, whose parameters are consistent with, for example, 
Moriya's self-consistent renormalization (SCR) theory of 
metallic spin fluctuations Jl3|-p^[ . 

We obtained these surprising conclusions and others to 
be discussed below through a neutron scattering study of 
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FIG. 1. The composition- pressure-temperature phase di- 
agram of the V2O3 system ^^6[||| . There exist four 
phases: the paramagnetic insulator (PI), the antiferromag- 
netic insulator (AFI), the paramagnetic metal (PM), and 
the spin-density- wave (SDW), a.k.a. antiferromagnetic metal 
(AFM). The Mott transition between the PM and PI phases 
is first order, as well as the transitions to the AFI phase. The 
first order nature of the latter phase transition becomes much 
less pronounced for Cr concentrations exceeding a;~3%. The 
transition to the SDW is second order. Except for the AFI 
phase which is monoclinic, all other phases have the trigonal 
corundum structure. 

magnetic correlations in V2O3 and its doped derivatives 
which constitute a famous 3-D Mott system |l9],fi] . Even 
though magnetic order occurs at low T in the metallic 
and the insulating phases of V2O3, the material is the 
only known transition metal oxide to display a Para- 
magnetic Metal (PM) to Paramagnetic Insulator (PI) 
transition (see Fig. [j]). The PM to Antiferromagnetic 
Insulator (AFI) transition in pure V2O3 is a spectac- 
ular first order phase transition in which the resistiv- 
ity abruptly increases by 8 orders of magnitude pO|pl[ , 
the lattice structure changes from trigonal to monoclinic 
^2j,^3|, and the staggered moment immediately reaches 
over 80% of the low temperature ordered moment p4[ . 
The mechanism for the phase transition remains contro- 
versial though it is clear that electron correlations are 
important because band structure calculations predict a 
metallic state for both the trigonal and monoclinic crystal 
structures. [p5|,p6[ One explanation for the transition is 
due to Slater }27|] who argued that the antiferromagnetic 
transition doubles the unit cell which in turn opens a gap 
at the Fermi energy. This mechanism however cannot ac- 
count for the PM-PI transition and offers no explanation 
for the existence of the PI phase. Another viewpoint is 
that the AFI phase is simply the magnetically ordered 
state of the PI and that the Mott mechanism works for 
the PM-AFI transition as it does for the PM-PI tran- 
sition. Following this approach, there have been vari- 
ous attempts to incorporate the antiferromagnetic phase 
in the framework of the Mott metal-insulator transition 



2 



p8| ^o|,^2|. However, none of the theories convincingly 
accounts for the first order nature of the PI-AFI tran- 
sition observed in V2O3. Another shortcoming of the 
theories is that they assume that low energy and static 
magnetic correlations are characterized by the same wave 
vector throughout the phase diagram, something that our 
experiments show is not the case for the V2O3 system. 

Magnetic order and spin waves below 25 meV in the 
AFI phase of V2O3 were previously measured by other 
workers [ pipHJ^l . Our experiments are the first neutron 
scattering study of magnetic correlations in the AFM, 
PM and PI phases of V2O3. Our most important con- 
clusions are the following: (1) the AFM has a small mo- 
ment incommensurate spin-density-wave (SDW) at low 
temperatures, which results from a Fermi surface insta- 
bility. (2) Dynamic spin correlations in the entire PM 
phase are controlled by the same Fermi surface instabil- 
ity. (3) Magnetic correlations in the metallic phases are 
those one might expect near a quantum critical point. 
They are quantitatively described by the SCR theory for 
itinerant antiferromagnetism in the small moment limit. 
(4) Anomalously short-range dynamic spin correlations 
in the PI are closely related to those in the PM but are 
different from the magnetic order in the AFI. (5) Phase 
transitions to the AFI, from either the PM or the PI, 
cause an abrupt switch of the magnetic wave vector, sig- 
naling a sudden change of the exchange constants in the 
spin Hamiltonian. All of (l)-(5) require more than a 
conventional Hcisenberg spin Hamiltonian for a satisfac- 
tory explanation. In metallic samples the geometry of 
the Fermi surface directly affects magnetic correlations, 
while in insulating samples there are low energy orbital 
degrees of freedom which are important for describing 
the unusual PI phase and the PI-AFI phase transition. 
Some of the above results were previously published in 
short reports |3^-|3S|]. 

The organization of the paper is as follows. Section 
II covers experimental details concerning sample prepa- 
ration and neutron scattering instrumentation, sections 
III, IV, V and VI describe magnetic correlations in the 
AFM, PM, AFI and PI phases respectively, and section 
VII concludes the paper with a discussion of the overall 
interpretation of the data. 

II. EXPERIMENTAL DETAILS 

V2O3 has the corundum structure |4(i[ ] (space group 
No. 165, R3c, refer to Fig. |) in all but the AFI phase, 
where small lattice distortions break the three fold rota- 
tion symmetry to yield a monoclinic structure. [ p3[ We 
use the conventional hexagonal unit cell |3^j2p[] with six 
V2O3 formula units per unit cell to index real and recip- 
rocal space. The reciprocal lattice parameters are a* = 
An/^/ia = 1.47(1)1" 1 and c* = 2n/c = 0.448(1)1^ 
in the metallic phases, a* = 1.46(1)A _1 and c* = 
0.449(2)A" 1 in the AFI phase, and a* = 1.45(1)A _1 and 




FIG. 2. The corundum lattice structure of V2O3. The solid 
symbols denote V ions and the hatched symbols O ions. 

c* = 0.451 (2) A" 1 in the PI phase. More accurate values 
of lattice parameters, as functions of temperature and 
doping, can be found in Ref. [ [^2|,|lJ^|]]. With this unit 
cell, the characteristic wave- vector of the AFI spin struc- 
ture is (1/2,1/2,0) while it is (0,0,1.7) for the SDW 
structure |33| (refer to Fig. ^|). Samples were oriented 
in either the (hhl) or the (hOl) zone to reach peaks in 
the neutron scattering associated with the two types of 
magnetic correlations. Selection rules for nuclear Bragg 
points in these two zones are (hhl): I — 3n and (hOl): 
h — I = 3ri, I = In. 

Single crystals of Ni^yO^ were grown using a skull 
mcltcr [Q. The stoichiometry of pure V2O3 and metal- 
deficient V2-J/03 was controlled to within Sy = 0.003 
by annealing sliced, as-grown V2O3 crystals in a suit- 
ably chosen CO-CO2 atmosphere Q for two weeks at 
1400°C. Single crystals of (Vi-zCr^Os were grown us- 
ing the Tri-arc technique ]2l] ]. The compositions were 
determined by atomic absorption spectroscopy and wet 
chemical titration technique. The mass of samples for 
elastic neutron scattering was typically 100 mg, while we 
used several grams for inelastic neutron scattering ex- 
periments. Large single crystals were sliced into plates of 
*~1.5 mm thickness as shown in Fig. |i](a) to allow oxygen 
diffusion to the bulk during the annealing process. Sev- 
eral single crystals of each composition were mutually 
aligned to increase the sensitivity for inelastic measure- 
ments. We report experiments for the seven compositions 
which are listed in Table Q. The samples have been well 
characterized inprevious studies using bulk measurement 
techniques |l]^|44|. 

The sample temperature was controlled using liquid 
He flow cryostats in the range from 1.4 K to 293 K, or 
using displex closed cycle cryostats from 10 K to 293 K. 
With a heating element, the displex cryostat could reach 
temperatures up to 600 K. 

Hydrostatic pressure up to 8 kbar was produced using 
a He gas high pressure apparatus [^0| at NIST, which can 
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(a) SDW (b) AFI 




FIG. 3. (a) Spin structure in the SDW state Q of 
V2-1/O3. Both spiral chiralities exist in our crystal. The stag- 
gered moment is 0.15/is per V. The four types of near neigh- 
bor spin pairs which are coupled by appreciable exchange in- 
teractions are indicated by A-D. (b) Spin structure in the AFI 
state [24J , which breaks the three- fold symmetry of the corun- 
dum structure. The staggered moment is 1.2/Lts per V. a-77 
indicate seven distinct exchange constants. The conventional 
hexagonal unit cell is used in this paper. The primitive vec- 
tors a and b of this cell are shown in (b), and the c vector, 
which is not shown, spans the stack from bottom to top. The 
primitive vectors of the primitive trigonal unit cell are shown 
in (a), and the primitive vectors a m , b m and c m of the prim- 
itive monoclinic magnetic unit cell of the AFI in (b). There 
are four spins (1-4) per primitive cell. 

TABLE I. Samples used in this study. For easy location of 
the samples in the phase diagram (Fig. yj) , x or y in molecular 
formula (Vi-^Cr^Oa or V2- B 03 is provided. The ambient 
transition temperature and physical states below and above 
the phase transition also are listed. 





y or x 


T N (K) 


K 


300 K 


Vl.96303 


0.037 


8.4(1) 


SDW 


PM 


V1.97303 


0.027 


8.5(1) 


SDW 


PM 


Vl.98303 


0.017 


9.5(1) 


SDW 


PM 


Vl.98503 


0.015 


55(1) 


AFI 


PM 


Vl.988 03 


0.012 


70(2) 


AFI 


PM 


V203 


0.0 


170(2) 


AFI 


PM 


Vl.944Cro.05603 


0.028 


180(2) 


AFI 


PI 



be loaded into an "orange" ILL cryostat with a temper- 
ature range from 1.4 K to room temperature. 



A. Neutron Scattering Instrumentation 

Neutron scattering measurements were performed on 
triple-axis spectrometers BT2, BT4 and BT9 at the NIST 
research reactor, on HB1 at the High Flux Isotope Reac- 
tor (HFIR) of ORNL, and on H4M and H7 at the High 
Flux Beam Reactor (HFBR) of BNL. For unpolarized 




FIG. 4. (a) Cross-cut pattern used to facilitate oxygen dif- 
fusion in single crystals during annealing, (b) Adjustable neu- 
tron aperture, made of Boron-polyethylene, to limit fast neu- 
trons within small angles of the main beam. Two pairs of 
opposite "jaws" can move into each other to adjust the cen- 
tral opening according to the sample size. 

neutron experiments, we used a Be (101) monochroma- 
tor on HB1, Cu (220) and pyrolytic graphite (PG) (002) 
monochromators on BT4 and PG (002) monochromators 
on all other instruments. PG (002) analyzers were used 
for all unpolarized triple- axis experiments. For polarized 
neutron measurements on BT2, crystals of Heusler alloy 
CU2M11AI [d(lll)=3.445A] were used as monochromator 
and analyzer. Fast neutrons were removed from the inci- 
dent beam by an in-pile sapphire filter on HB1, and 
PG filters ]52[ | were used to remove high-order neutron 
contaminations where appropriate. The horizontal colli- 
mations were controlled with Soller slits and are specified 
in the figures containing the experimental data. The fast 
neutron background at BT4 in the small scattering an- 
gle limit was reduced by a 30 cm long boron-polyethylene 
aperture before the sample and a similar 15 cm long aper- 
ture after the sample. The opening of the aperture was 
adjusted to match the sample size (Fig. ||(b)). These 
apertures suppress fast neutron background in the small 
angle limit which is important for measuring high energy 
inelastic magnetic scattering. 

In most of our scans, the background is approximately 
constant. For some scans, however, the scattering angle 
falls below w 7° where background increases rapidly with 
decreasing angle. For example, a constant-ftw = 28meV 
scan measured with H4M is shown in Fig. || with solid 
circles. The scattering angle at (100) is 5.4°. As a result, 
the background (open symbols), measured with the ana- 
lyzer turned 10° from the reflection condition, increases 
near (100). The solid line represents data with this extra 
background contribution subtracted |53| . 

B. Neutron Scattering Cross Section 

The magnetic neutron scattering cross-section for mo- 
mentum transfer hq = h(k — k') and energy transfer Tiuj 
is given by [|| 
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V i era 3 at 150K, ft« = 28meV 




-2-101234 
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FIG. 5. A constant energy scan (solid circles) for which the 
scattering angle near (100) is close to the incident neutron 
beam. The squares denote the background measured with 
the analyzer rotated 10° from the reflection condition. See 
text for details. 



The intensity in a neutron scattering experiment is 
measured as the ratio of neutron counts in a high effi- 
ciency 3 He detector to the neutron counts in a low ef- 
ficiency monitor placed in the incident beam before the 
sample. To cover a wide energy range of spin excitations, 
several spectrometer configurations were used. The rela- 
tive sensitivity of different configurations was determined 
by comparing the integrated intensities of identical scans 
probing magnetic or incoherent nuclear scattering. Nor- 
malizing to inelastic scattering from acoustic phonons 
p6| , p7t yielded absolute measurements of <S(q, u) [Eq. (||)] 
in units of /i^/meV. Units of mbarn/meV are also used 
in this paper and when we do so we are quoting numbers 
for the following normalized intensity 

/(q, W ) = 2(^) 2 [/(q)] 2 5(q,c), (6) 

which is more directly related to the raw data because it 
contains the magnetic form factor, /(q), as a factor. 



a/3 



(1) 



where N is the number of hexagonal unit cells, 
( 7 r /2) 2 = 72.65 x 10" 3 barn/u|, /(q) is the magnetic 
form factor for the V 3+ ion |55[ , exp[— 2M / (q)] is the 
Debye- Waller factor, and the dynamical structure factor 
is given by 

^ (qiW) = MiV /"^M-q-CR-R')] 



RR' ' 



x(S£(i)S&,(0)). 



(2) 



g is the Lande factor, S^(t) is the ath Cartesian spin 
component of the V ion at position R at time t, and (. . .) 
denotes a thermal average. When spin space anisotropy 
can be neglected, as in the paramagnetic phase of a 
Heisenberg system, expression (Q) simplifies to 
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^=4 (f0 2 ^ 2e ~ 2W(q) ^ aQ ^ 



(3) 



We will drop the superscripts in the remaining text for 
brevity. The structure factor is related to the imaginary 
part of the dynamic spin susceptibility by the fluctuation- 
dissipation theorem 
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^ Q—Hio/kt 



rV(q,w). 



(4) 



The real part of the generalized spin susceptibility is 
given by the Kramers-Kronig relation, 



x'(q,w) - x'(q,°°) 



du) 



,x"(q,w') 



(5) 



C. Resolution Effects 

The measured intensity in a neutron scattering experi- 
ment is proportional to the convolution of the scattering 
cross-section with an instrumental resolution function. 
The resolution function can be approximated by a Gaus- 
sian function in 4-dimensional q-w space 1 58 5^] . Inelastic 
neutron scattering experiments are usually flux-limited. 
Therefore the choice of an experimental configuration for 
probing a specific scattering cross section always involves 
a trade-off between resolution and sensitivity. The fol- 
lowing two limits are of importance in this study. 



1. Infinite life time excitations 

Examples of this type of excitations, such as phonons 
and un-damped antiferromagnetic spin waves, are de- 
scribed by a scattering cross section of the form 

5(q,w) ~ 5(uj - c|q|). 

If the curvature of the dispersion surface is negligible over 
the volume of the resolution function then the dispersion 
surface can be approximated by a flat plane and the re- 
sulting profile of a scan through the dispersion surface 
is approximately Gaussian with a width defined by the 
resolution. With coarser resolution, the profile becomes 
asymmetric and the peak position may shift from the 
position defined by the dispersion relation. Refer to the 
dotted lines in Fig. nfl for a few examples of such effects. 



2. Weakly q and w dependent scattering cross section 

Features in the scattering cross-section which are 
broader in q-u> space than the experimental resolution 
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are reproduced directly or with slight additional broad- 
ening in the measured q and lo dependent intensity. This 
is the case for our measurements in the AFM, PM and 
PI phases. For example if 6>(q, u>) takes the form 



V 1.983°3 at !- 4K 

14.flmeV, BQ-Heusler-fJ0-60-Heusler-6O f PG filter 



5(q, 



exp(-|q-q | 2 /crg), 



a constant energy scan will yield a Gaussian peak with 
width (ctq -t-cfj) 1 / 2 , where or is the resolution width. To 
appreciate the resolution broadening, let us put in num- 
bers: When the resolution width is 1/3 of the intrinsic 
width, the measured width increases by less than 10% 
over the intrinsic width. When the resolution width is 
1/5 of the intrinsic width, the measured width is only 
2% larger than the intrinsic width. In this latter case the 
measured intensity is then very close to being directly 
proportional to 5(q, to). 



III. SDW IN METALLIC V 2 3 

The insulating antiferromagnetic state (AFT) of V2O3 
can be suppressed by a pressure of ~20 kbar |6Q|JXs| ] - The 
pressure-stabilized metallic state is paramagnetic down 
to at least 0.35 K. The AFI state can also be 

suppressed at ambient pressure by vanadium deficiency 
[]62fl or titanium doping |19| . The critical concentra- 
0.1 for Va-^Ti^Oa 



19| and y c ~ 0.015 

The existence of antiferromagnetic or- 



tions are x c 
for V 2 _,0 3 |5[ 
der at low temperatures for the doping induced metal 
was first established through the 57 Fe Mossbauer effect 
|p3] p5[ . The magnetic structure was only recently de- 
termined through single crystal neutron diffraction [ p3[ , 
and it is depicted in Fig. [}] (a). This magnetic order can 
be described as puckered antiferromagnetic honeycomb 
spin layers whose spin directions form a helix along the 
c axis. The spin directions determined from our neu- 
tron diffraction refinement, which lie in the basal plane, 
are consistent with previous susceptibility anisotropy and 
magneto-torque measurements J64| , p6| . The pitch of the 
helix yields an incommensurate magnetic structure, with 
magnetic wave vector 1.7c*. Both the staggered moment 
and the wave vector depend weakly on doping, y. |3^] 
This magnetic order in metallic V2- y 03 bears little re- 
semblance to the magnetic order in the insulating phase 
of V2O3 which is characterized by a magnetic wave vector 
(1/2,1/2,0) H [refer to Fig. 0(b)]. 

We have previously shown ]33[ ] that the antiferromag- 
netic order in the metallic state can not be described by 
a localized spin model, since the dynamic magnetic spec- 
tral weight greatly exceeds the square of the staggered 
moment and the bandwidth for magnetic excitations ex- 
ceeds ksTjv by more than one order of magnitude. The 
proper concept is probably that of a spin density wave 
p7| p9[ where magnetic order results from a Fermi sur- 
face nesting instability. This interpretation is supported 
by recent RPA calculation of the q-dcpendcnt suscep- 
tibility using realistic band structure [M. According to 
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FIG. 6. Triple-axis polarized neutron scans through the 
(1,0,0.3) magnetic Bragg peak with neutron spins polarized 
along the momentum transfer (HF) in (a), and along the ver- 
tical direction (VF) in (b). Measurements with the neutron 
spin flipper on and off are presented. The contribution to 
scattering from spin components in the basal plane is shown 
in (c) and the contribution from components polarized along 
the c axis is shown in (d). See text for details. 



neutron, specific heat, and transport measurements, only 
a small fraction of a large Fermi surface |33| is involved 
in the SDW of metallic V 2 - y 3 . As seen at the SDW 
transition for Cr metal Jfl[ , the resistivity for metal- 
lic V2_ a C>3 increases below the Neel temperature as the 
Fermi surface is partially gapped. Q This is in contrast 
to localized spin systems where the resistivity decreases 
below Tn because of the decrease in spin disorder scat- 
tering. The discovery of the transverse incommensurate 
SDW resolved the longstanding mystery Jt^] about the 
magnetic ground state of doped metallic V2O3. 



A. Polarized neutron scattering 

To further characterize the incommensurate magnetic 
structure, polarized neutron measurements were per- 
formed. A V1.983O3 sample was oriented with its (hOl) 
zone in the horizontal scattering plane. A small magnetic 
field of 7.6 G was applied along either the scattering wave 
vector (HF) or along the vertical direction (VF) to guide 
the neutron spin. A neutron spin flipper was inserted in 
the neutron path. In the HF case when the spin flipper 
is turned on, only magnetic scattering contributes to the 
Bragg peak |73J . Coherent nuclear scattering is non-spin- 
flip and contributes when the flipper is off. Therefore, 
Fig. ||(a) shows that the (1,0,0.3) Bragg peak is indeed 
magnetic. The finite intensity in the flipper off case is 
due to incomplete polarization of the neutron beam and 
is consistent with a flipping ratio of /y=14, measured at 
nuclear Bragg peaks. 
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Since 0.3c*/a* « 0.09, (1,0,0.3) is essentially parallel 
to the a* axis. The four partial intensities at (1,0,0.3) for 
the HF and VF cases with the spin flipper on or off are 
therefore Iral 



jon 
1 HF 

joff _ 



Ic, 



r& = {h + h)/fH, 

Ivf = Ic + h/fv, 

Iy" =h + Ic/fv, 



(7a) 
(7b) 
(7c) 
(7d) 



where L 



6(c) 



*6(c) 



is the contribution from the spin 
components in the b = a* x c* (or c) direction, and 
fmv) — 14 is the flipping ratio in the HF (VF) case. 
Therefore, from the data shown in Fig. || (a) and (b), 
and I c can be extracted by solving the over-determined 
equations (|7a])-([7d|). For example, 

I h = I° v f l - I c /f v 



Ic. = 



joff 
1 VF ' 
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I VF 



joff 
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In the last step in (g), a small quantity I c /fv is ne- 
glected. Alternatively, 
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Fig. ^ (c) and (d) show the basal plane component and 
the vertical component I c derived from the data in Fig. ^ 
(a) and (b) using these expressions. The solid circles 
were determined using Eq. (ph and the open circles using 
Eq. (j^). The two data sets are consistent within exper- 
imental uncertainty and show that the incommensurate 
peak is associated with basal plane spin correlations as 
was previously conjectured on the basis of unpolarized 
neutron diffraction. The data also reveal a second in- 
commensurate peak with mixed polarization. This peak 
has only been observed in this sample which is very close 
to the critical concentration, y c 0.015, for the ambient 
pressure metallic phase J|3| . 

When polarized neutrons are scattered by a right- 
handed transverse spiral in the HF configuration, the 
partial cross-sections for magnetic Bragg peaks with wave 
vectors along the spiral axis are 1 73 



do-(q) 
dQ 



H — »- 



oc^e 



i(q+q m ).(R-R') 



RR' 



^<5(q + q m - r) 



(10) 



14.8meV, 60-Heusler-60-samp!e-60-open, PG filter 



1000 



C 

«J 500 



i — 1 — i — r ~ 

HF: • on 

o off 



i L 



500 



I(1.4K)-I(15K) HF: . on 
o off 




7.4 7.6 7.8 B.O 



7.4 7.6 7.8 8.0 



(00 /!) 



FIG. 7. (left) Double-axis polarized neutron scans at 
(0,0,7.7) with the neutron spin polarized along the c axis (HF) 
at 1.4K. (right) Identical scans at 15K were subtracted as a 
measure of non-magnetic background. The magnetic peaks 
are independent of the incident neutron spin state which 
shows that the SDW phase has no macroscopic handedness. 



for the neutron spin up to spin down channel and 



<Mq) 



RR' 

^(S(q-q TO -r) 



(11) 



for the neutron spin down to spin up channel. In these 
expressions the nuclear reciprocal lattice vector r and 
the spiral wave vector, q m , are both parallel to the spi- 
ral axis. Note that only one of the two magnetic satel- 
lite Bragg peaks appears in each channel. For a left- 
handed spiral, the signs of q m in these expressions re- 
verse. Therefore if the SDW in V2- y C>3 is a spiral with 
a macroscopic handedness, we should find that magnetic 
Bragg peaks with Q = (00/) come and go with a 180° 
rotation of the incident neutron spin state. 

Fig. ^ shows a double-axis polarized neutron scan 
through the (0,0,7.7) magnetic Bragg peak of the SDW. 
The left panel shows raw data taken at 1.4 K, while in 
the right panel, a non-magnetic background taken above 
the Neel temperature has been subtracted. Clearly the 
magnetic Bragg intensity does not change upon flipping 
the incident neutron spin state. Hence we conclude that 
the SDW phase does not have an intrinsic macroscopic 
handedness. This means that if the ordered phase is a 
spiral, our sample contains equal volume fractions of left 
and right handed spirals. Another possibility is that we 
have an amplitude modulated SDW which is a micro- 
scopic superposition of the left and right-handed spirals. 



B. SDW, PM and AFI phase boundary at T=0 

After confirming the magnetic structure in the AFM 
phase with polarized neutrons, let us now try to delimit 
the phase boundary of this SDW ground state in the 
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FIG. 8. Phase transitions at high pressure for V1.988O3. 

(a) Suppression of the AFI state by pressure, as measured by 
monitoring the (1/2,1/2,0) magnetic Bragg peak upon cool- 
ing. Insert shows scans through (1/2,1/2,0) at 3 kbar above 
(open circles) and below (solid circles) the AFI transition. 

(b) Elastic scans at 1.4 K and 5.0 kbar, covering a region in 
the (hhl) zone slightly larger than a quadrant of the Brillouin 
zone. Besides two spurious peaks which were found to be tem- 
perature independent, there are only two resolution-limited 
nuclear Bragg peaks, (110) and (113). [No data were taken at 
the two corners near (1/2,1/2,3) and (113).] (c) Elastic neu- 
tron scans at 1.4 K through (1,0,2.3) under pressures which 
stablizes the AFI (solid circles) and the SDW (open circles) 
ground states respectively, (d) The pressure-temperature 
phase diagram for V1.988O3. The hatched region shows the 
thermal hysteresis at this first order transition. 



pressure-composition plane for V2- a C>3. Two of the sam- 
ples we studied with high pressure neutron diffraction 
(y = 0.012 and 0.015) have the AFI ground state at am- 
bient pressure and display the metallic state only under 
pressure. Fig. || (a) shows the (1/2,1/2,0) Bragg peak 
which is the order parameter of the AFI, as a function of 
temperature at various pressures for a V1.988O3 sample. 
Under high pressure, the transition to the AFI remains 
strongly first-order, as shown here by the sudden onset of 
antiferromagnetic order with a nearly saturated magnetic 
moment. 

The AFI state is completely suppressed above 4.1 kbar 
for V1.988O3. A survey at 5 kbar and 1.4 K in a quadrant 
of the (hhl) Brillouin zone reveals no magnetic peaks 
associated with the AFI phase [refer to Fig. || (b)] . On the 
other hand the Bragg peaks associated with SDW order 
are absent for P=4.0 kbar but appear under a pressure of 
4.6 kbar [refer to Fig. || (c)]. The temperature-pressure 
phase diagram for V1.988O3 is shown in Fig. |§| (d). The 
SDW transition is second order, |33||36| while the first 
order AFI transition shows large thermal hysteresis. 

Fig. |^ summarizes our current knowledge about ground 
states in ¥2-3,03. For y = 0, after the AFI is suppressed, 
the metallic state is paramagnetic. For y = 0.012 and 
0.015, we found that after the AFI is suppressed, the 
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FIG. 9. The composition-pressure pha se diagram at T = 
for V2- H 03. The squares are from [ and circles from 

the current neutron mesurement. Open circles are extrapo- 
lated from data below 8 kbar. For y ~ 0, the paramagnetic 
metallic (PM) state is stable at high pressures. However, the 
extent of the PM state beyond the hatched region and down 
to temperatures below 0.1 K is not clear. It has not been pos- 
sible to produce single phase samples of V2- y 03 for y > 0.05 
because of the formation of V3O5 



metallic state is a SDW. Up to 8 kbar, the upper limit of 
our pressure cell, the SDW remains the stable low tem- 
perature state for all the samples (y = 0.012, 0.015 and 
0.027) which we have studied with high pressure neutron 
diffraction |33|]36] ]. By linear extrapolation of the pres- 
sure dependent Neel temperatures, we derived estimates 
for the critical pressures of the SDW phase and these are 
indicated by open circles in Fig. ^. 

Anderson has pointed out the possibility of a super- 
conducting instability in the Brinkman-Rice liquid [fT4[ . 
For metallic V2~yOs, SDW order appears to be the dom- 
inant instability. However, it should be noted that the 
possibility of a superconducting state with T c < 0.1 K at 
high pressures has not been excluded yet WB. 



IV. SPIN EXCITATIONS AND QUANTUM 
CRITICAL BEHAVIOR IN METALLIC PHASES 

A. Stoner electron-hole pair continuum at T < Tjv 

Fig. |l^ shows constant energy transfer, hui, scans 
spanning a whole Brillouin zone in the AFM phase for 
V1.973O3 (Tw = 8.5(1) K). The incommensurate Bragg 
points, (1,0,0.3) and (1,0,2.3), covered by these scans 
correspond to the long range ordered SDW. The solid 
circles represent a flat background measured by turning 
the PG(002) analyzer 10° away from the reflection con- 
dition. Around each magnetic Bragg point, (1,0,0.3) or 
(1,0,2.3), there is only one broad peak. This is not due 
to poor resolution merging two counter-propagating spin 
wave modes: the dotted lines in the figure represent the 
calculated spectrometer response to spin waves with a 
velocity c = 130 mcVA, which is much larger than the 
inverse slope of the peak width versus energy transfer 
curve (see Fig. O). Our resolution is clearly adequate 
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FIG. 10. Constant energy, Tiuj, scans at 1.5 K for V1.973O3 
across a Brillouin zone. (102) and (104) are two nuclear 
Bragg points, and (1,0,0.3) and (1,0,2.3) are two incommensu- 
rate magnetic satellites (marked by crosses). The spectrom- 
eter configurations at BT2 were Ef = 13.7 meV with hor- 
izontal collimations 60'-40'-40'-40' for (a)-(c), and Ei = 35 
meV with 60'-40'-60'-60' for (d). Vertical collimations were 
300'-170'-230'-970' in all four cases. The filled circles are an- 
alyzer turned backgrounds. The dotted lines represent the 
four dimensional convolutions of the instrumental resolution 
with an isotropic spin wave dispersion relation with a velocity 
c = 130 meVA. 



to rule out a conventional spin wave response in this en- 
ergy range. Therefore the broad intensity profile which 
we observed is intrinsic, and in the energy range probed, 
magnetic excitations consist of a single broad lobe at each 
magnetic Bragg point. 

For insulating magnets, propagating spin waves exist 
over the whole Brillouin zone. For itinerant magnets, 
however, propagating spin waves only exist in a lim- 
ited range near the magnetic Bragg peaks. Beyond this 
range they merge into the Stoner particle-hole contin- 
uum. According to weak coupling theory, |6^] the energy 
gap which a SDW creates inside the particle-hole contin- 
uum is A = 1.76afcsTjv, where a = (yi + V2)/2{yiV%) l l 2 
and V\ and v 2 are Fermi velocities of the two nesting 
bands. Inside this gap, spin waves propagate with a ve- 



locity @7§|73] c = (viv 2 /3) 1/2 . For metallic V 2 ~ y 3 at 



T=0, spin waves are expected only below ~1.76fc.gT/v = 
1.4 meV, which is below the energies which we have 
probed in this experiment. The broad lobes measured 
in Fig. [l^ thus represent magnetic excitations inside the 
Stoner continuum. 

Other signs that we are probing magnetic fluctuations 
in the electron-hole continuum, as opposed to transverse 
oscillations of the order parameter, are (1) The spectral 
weight from 1.5 meV to 18 meV in Fig. [10| corresponds to 
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FIG. 11. Resolution-corrected peak widths (FWHM) of 
constant energy scans as a function of energy transfer for 
V1.973O3 in the long range ordered SDW state at 1.5 K. 



a fluctuating moment of 0.32/i# per V, already twice that 
of the static staggered moment, and (2) The resolu- 
tion corrected full-widths-at-half-maximum (FWHM) of 
constant energy scans (see Fig. [ll]) extrapolate to a finite 
value at to — 0, even though there are true, resolution- 
limited, Bragg peaks at co = 0. 

For Cr metal another 3-D SDW (T N = 310 K), 

one would expect at T=0 that spin wave modes merge 
into a single broad lobe above ~1.76A;b2jv=47 meV. 
However, due to the much lar ger sp in wave velocity in 
Cr, estimated at ^1000 mcVA |8C|j8l]] , it has not yet been 
possible in this material to distinguish between resolution 
merged spin waves and single lobe excitations. We were 
aided in the case of metallic V2O3 by strong correlations. 
The enhanced effective mass reduces the Fermi velocity 
to a value which is ~ 1/15 of that for Cr. |53| This makes 
it feasible to unambiguously observe these unique single- 
lobe magnetic excitations in the SDW state of V2-j,03. 



B. Doping effects in the metallic phases 

We have previously shown |3^] that both the stag- 
gered moment and the incommensurate wave vector of 
the static SDW order depend only weakly on doping. For 
a change in metal-deficiency Sy = 0.03 which is twice the 
critical y c for the AFI phase, the change in the magnetic 
wave vector (1.7c*) is only ~0.02c*. Likewise, dynamic 
magnetic correlations in the PM phase are insensitive 
to doping. As shown in Fig. |l2|, even for pure V2O3 
which has different low temperature AFI state, the mag- 
netic response at 200 K is indistinguishable from that 
of V1.973O3. In particular we do observe peaks at the 
SDW wave vectors. This demonstrates that magnetism 
throughout the metallic phase is controlled by the same 
Fermi surface instability and indicates that small mo- 
ment incommensurate antiferromagnetism in V2-yC>3 is 
not due to impurity effects but is an intrinsic property of 
metallic V2O3. 
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FIG. 12. Magnetic responses at 200 K from V1.973O3 (open 
circles) and V2O3 (solid circles) in the PM phase. Although 
these two samples have different low temperature states, the 
SDW for the former and the AFI for the latter, there is no 
difference in the dynamic correlations in the PM phase. The 
experimental configuration was E/ = 13.7 meV with hori- 
zontal collimations 60'-40'-40'-60' for V1.973O3 at BT2. The 
horizontal bars indicate the FWHM of the projection of the 
resolution function on the scan direction. The experimental 
configuration for V2O3 were E/ = 13.7 meV with horizon- 
tal collimations 40'-60'-80'-80' for Hlu = 6 and 9 meV, and 
Ei = 35meV with 40'-40'-80'-80' for 12 meV at BT4 with PG 
monochromator and analyzer. The solid line is a least-squares 
fit to Eq. and @. 



C. Antiferromagnetic quantum critical behavior in 
metal 

The SDW in metallic V2-j,03 is a small moment, low 
temperature antiferromagnetic ordering. The staggered 
moment in the magnetic ground state (0.15/is) is only 
a small fraction of the effective paramagnetic moment 
(2.37-2.69/is) which appears in the high temperature 
Curie-Weiss susceptibility j8^^|j(|. In other words, the 
magnetic system is very close to a quantum critical point 
[ pjj p7{ . We have previously shown |Q that the sin- 
gle lobe spin excitation spectrum in metallic V1.973O3 is 
isotropic in space, and the spectrum x"(q, from 1.4K 
inside the ordered phase to temperatures above 20Tjy, 
has been determined in absolute units [Q. Here we dis- 
cuss our results in connection with the quantum critical 
behavior. 



Fig. 13(a) shows in a log- log scale the inverse correla- 
tion length, k, measured in Rcf. [ |$(|] as a function of the 
reduced temperature, t = (T — Tn)/Tn- The solid line 
is k oc t v with v = 0.55(4) from a least squares fit. This 
value of v suggests a small positive r\ exponent according 
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FIG. 13. (a) The inverse correlation length as a function of 
the reduced temperature for V1.973O3 in the PM phase. Both 
axes use the logarithmic scale. The solid line is the best fit of 
K oc t" to the data with v = 0.55(4), and the dotted line is a 
fit with fixed v = 1/2. (b) The characteristic spin fluctuation 
energy, hui p , where ^"(Q,^) is maximized, as a function of 
k. The solid line is the best fit of u p cx k z to the data with 
z — 1.9(1), and the dotted line is a fit with fixed z — 2. (c) 
hu p is roughly proportional to the thermal energy ksT (the 
dotted line). The solid line is ui p oc T — T/v. 



to the scaling relation p9| 77 = 2 — 7/1/. However, the 
classical values of v = 1/2, 77 = are also consistent with 
our data (refer to the dotted line). 

The spin fluctuation spectrum %"(q, w) reaches its 



and q = Q which is 
The dominant 



maximum value, x'maxi at w = 
an antiferromagnetic Bragg point 36 
part of the spectral weight, therefore, falls in the hydro- 
dynamic region |m|, |q — Q| -C k, where fiui p — jak z . 
In Fig. |l3](b), we find the dynamic exponent z = 1.9(1) 
(the solid line), which is in agreement with the theoreti- 
cal prediction z — 2 for itinerant antiferromagnets in the 
quantum critical region [ p4| , ^7| . The constant prefactor 
is 7a = 87(4) meVl 2 . 

Combining the scaling relations from Fig. |l3|(a) and 
(b), we have huj p ~ t zv = CksiT — Tv) [refer to the solid 
line in Fig. |l3](c)]. Given the small value of Tn and the 
fact that C ~ 1, the characteristic spin fluctuation energy 
fiLOp is close to fcsT [refer to the dotted line in Fig. 13(c)!. 
This feature is shared with nearly metallic cuprates |9l| ] 
and some heavy fermion systems 1 92 . Antiferromagnetic 
spin fluctuations for systems close to a quantum critical 
point, therefore, might provide a microscopic mechanism 
for marginal Fermi liquid phenomenology |)3j . 

The peak intensity of some constant energy scans is 
shown in Fig. |lj(a) as function of energy. They are ex- 
pressed in absolute units and represent x"(Q,w) at the 
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FIG. 14. (a) Normalized peak intensity of constant energy 
scans at 1.5K (triangles), 50K (diamonds), 130K (squares) 
and 200K (circles) for V1.973O3. In the measured energy 
range, x"(Qi w ) changes from a decreasing function of Tilu 
at 1.5K to an increasing function at 200K. (b) Data in 
(a) are scaled onto a universal curve using measured 
Xmax = Xq/2 and haj p = 7 a k 2 ■ The open circles are from 
constant-?ic 1 ;=1.5meV scans at 9K, 19K, 27K and 82 K; the 
open diamonds are from a constant-Q=(l,0,2.3) scan at 30 K. 
Constant-Q scan for V2O3 at 200 K is shown with the open 
squares. The solid line is z = 2y/(y 2 + 1) [refer to Eq. (113) 
and (HI)]. 



magnetic zone center. At 1.5K, x"(Q,cj) is a decreas- 
ing function of u> in the measured energy range and it 
changes to an increasing function at elevated tempera- 
tures. They were scaled onto a single universal curve in 
Fig. 0(b) by dividing x"(Q,w) by Xmax an d Hu by the 
peak energy %uj p . The apparently different behaviors in 
a finite energy window in Fig. |l4|(a), therefore, reflect 
different segments of a universal behavior. 

For a general wave vector, q, away from an antifer- 
romagnetic Bragg point, Q, we have shown [|36| previ- 
ously that our neutron scattering data measured at var- 
ious temperatures can be described by 



X "(Q + q, W ) 



2* 



max 



K K 



where the scaling function is 

$(x,y)=y/[y 2 + (l + x 2 ) 2 ] 



(12) 



(13) 



The solid line in Fig. [TJ(b) is 2 $(0,y). The maximum 
Xmax i s related to static staggered susceptibility fl36| . 
Combining Kramers-Kronig relation and Eq. (|lj) gives 



/oo 
du>x"(Q,u)/u 
-OO 

= X ' max 2TT- 1 / ^>{Q,hLo/ lA K 2 )dLo/u: 



$(0,y)dy/y 



X, 



x constant. 



Notice that when Tjv — > 0, = Cfc^T, Xmaa; 
aiT~ 7 and k = T" /a 2 , hence Eq. (p"2|) becomes 



ai 



x "(Q + q,c) = ^<I> 



a 2 q Tiw 



T" ' Ck B T 



(14) 

XQ = 

(15) 



where a\ and a 2 are non-universal constants. p4fl Using 
zzv fa 1 from our measurement and the scaling relation 
7 = (2 — n)i>, and absorbing the constant C in the scaling 
function, 



X "(Q + q,^) 



ai aj^*,^). (16) 



j>{2-r,)/z 



T 1 /^ k R T 



This expression, which is supported by our experimental 
data, is identical to the scaling form proposed by Sachdev 
and Ye in the quantum critical region [p6| . 



D. Theory of itinerant magnetism 

Itinerant magnetic systems have been studied within 
the random-phase approximation (RPA), which provides 
a reasonable description at T=0 for the ordered moment 
and for spin excitations in the spin wave modes and in 
the Stoner continuum. But it does not relate the low- 
temperature excitation parameters to finite-temperature 
properties such as the ordering temperature and the high 
temperature Curie- Weiss susceptibility. 

The imaginary part of the RPA generalized suscepti- 
bility for an itinerant antiferromagnet, as was done to 
describe magnetic fluctuations in chromium |7^,^], can 
be written as follows: 



X "(Q + q,^) 



Xo/«o) 2 hu}T(q,K) 



(17) 



l + (q/ K ) 2 (M 2 +r(q, K ) 2 ' 

where the wave vector dependent relaxation rate, 

r(q, K )= 7 ,4( K 2 + q 2 ), (18) 

corresponds to the simplest functional form consistent 
with the dynamical scaling hypothesis with a dynamic ex- 
ponent 2 = 2. As a comparison, for the insulating Heisen- 
berg antiferromagnet, |H)|,[)7| z = 3/2. The uniform non- 
interacting susceptibility xoj the length scale £oj and "/a 
are temperature-independent parameters, and the inverse 
correlation length k approaches zero when the Neel tem- 
perature is approached. This functional form and deriva- 
tions thereupon 198] offer a good description of critical 
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spin fluctuations for Cr and its alloys |8^,Q . They have 
also been successfully used to describe spin dynamics 
as measured by neutron scattering for Lai .sfiSrn.u CuCh 



V 3 3 in PM and AFI phases 



above the superconducting transition temperature [ 100 1 . 
A modified form was used to model nuclear magnetic 
resonance (NMR) dat a in othe r high temperature super- 
conducting cuprates [ 101 ~ 103| . However, the RPA pa- 
rameters 70, Xo, and £o in Eq. (17) are free parameters 
which are not connected to other physically measurable 
quantities. 

The SCR theory, |l4|-|l8| a mode-mode coupling the- 
ory, improves upon the RPA theory by considering self- 
consistently the effects of spin fluctuations on the thermal 
equilibrium state. It is solvable in the small moment limit 
with a few experimentally measurable parameters. The 
ferro mag netic version of the SCR theory was discussed in 
[ |]l^-|l7|l and the antiferromagnetic version in [ |l^ , |95[] 1. 

The generalized susceptibility for antiferromagnets in 
the SCR theory |Q has the same functional form as in 
the RPA theory: 



X "(Q + q,^) 



XQ 



hujT(q, k) 



l + (q//c) 2 (M 2 + r(q,K) 2 ' 



(19) 



where r(q, k) is defined as in Eq. (|18|). The RPA pa- 
rameter Xo/( K £o) 2 is substituted now by the static stag- 
gered susceptibility xq > m consistence with the Kramers- 
Kronig relation. Major difference from the RPA theory, 
however, is that the SCR theory can predict other physi- 
cal quantities using parameters from spin excitation spec- 
trum |pq| . For example, the non- universal parameters 
p4j a\ and 02 defined above Eq. ( |i~5| ) can be used to 
quantitatively relate Tjy to the staggered moment Mq 
(refer to Fig. 4 of Ref. [ 

Fig. [[5] shows an example of a simultaneous fit of 
constant energy scans to Eq. (|TJ) and ( |l8| ) for stoichio- 
metric V2O3 at 200K. Solid circles represent the magnetic 
response in the paramagnetic metallic phase, which peaks 
around the SDW wave vectors. The resolution widths 
are much smaller than the measured widths, therefore 
resolution broadening can be ignored in this case. The 
solid line through these 200K data is the best-fit func- 
tional form for |/(q)| 2 <S(q, oS) calculated using copies of 
Eq. (19) summed over the Bragg points Q = (1, 0, 4— A) 
and (1,0,2 + A) with A = 1.95(6)c*, « = 1.32(5)c*, 
7a = 18(2) meV/c* 2 and xq = Q.071(7)/4/meV per 
vanadium. These values are within the error-bars of those 
for V1.973O3, H^] as should be expected given the close 
similarity of the raw data for these two samples which 
are shown in Fig. |TJ[ 

It is customary to use the Lorentzian relaxation en- 
ergy T(q, k) of ( |i~7| ) or ( |j~9| ) as the energy scale of spin 
fluctuations. For an itinerant antiferromagnet close to a 
quantum critical point, such as metallic V2-y03, this 
reduces to jak 2 at q = To demonstrate the 

prominence of this energy, we plot in Fig. |lq(a) the in- 



tensity contour of x"(Q + W )/XQ = ^(l/K.hoj/^A n 2 ) 
for small moment itinerant antiferromagnets [refer to Eq. 
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FIG. 15. Constant energy scans at 200 K in the PM phase 
(solid circles) and at 160 K in the AFI phase (open circles) 
for stoichiometric V2O3. Magnetic correlations near the SDW 
wave vectors, represented by the two peaks, disappear when 
the material enters the AFI phase. The experimental config- 
urations at BT4 were E/ = 13.7 meV with horizontal colli- 
mations 40'-60'-80'-80' for Hlu = 6 and 9 meV, E; = 35meV 
with 40'-40'-80'-80' for 12 and 18 meV, and E f = 30meV with 
60'-60'-60'-60' for 25 meV. We used PG (002) as monochro- 
mator and analyzer except for the last configuration where 
a Cu (220) monochromator was used to improve the energy 
resolution. The horizontal bars indicate the FWHM of the 
projection of the resolution function on the scan direction. 
The solid line is a least-squares fit to Eq. ( |l9| ) and ( |l^ ) and 
absolute units are given on the right scales. 



(|i~9|), ( [l8| ) and (E3)]. The spectrum is unique in that the 
electron-hole damping has reduced the excitation spec- 
trum to a single lobe at the magnetic zone center, as we 
observed in our neutron scattering experiment for metal- 
lic V2- y 03 ]33| , |37| ] . The peak intensity occurs at q = 
and hu) = ta/? , and 7^4 n 2 ~ t zv clearly is the character- 
istic energy scale. In this case, zv w 1, as we found exper- 
imentally, is essential for achieving the quantum critical 
scaling of (|l6|). 

This is quite different from the well studied case of 
small moment itinerant ferromagnets. The x /; (q, ^) pre 



dieted by the SCR theory for ferromagnets [104,105] has 
the same Lorentzian form as (19), but the relaxation en- 
ergy is given by 



T(q, k) = 7F q(K 2 +q 2 ). 



(20) 



Like the antiferromagnetic SCR theory, the ferromag- 
netic SCR theory not only correctly describes the spin 
fluctuation spectrum, but also quantitatively relates the 
spectrum parameters to various other physical quantities, 
as demon strated i n the 3-D s mall-m oment ferromagnets 
MnSi and Ni 3 Al pit . Recasting Eq. © 

and (EQ) in a dimensionless form, 
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FIG. 16. Intensity contour of the scaled magnetic response 
function x"(Q + Qi^/xq as function of reduced wave vec- 
tor and reduced energy for (a) itinerant antiferromagnets and 
for (b) itinerant ferromagnets near a quantum critical point. 
The contour levels are 0.1, 0.2, 0.3 and 0.4. In the antifer- 
romagnetic case (a), the peak intensity is 0.5 at q=0 and 
Tuu/'yAii 2 ~i- In the ferromagnetic case (b), the intensity is 
singular at q=0 and uj—0. Both k and XQ are functions of 
temperature. 7a and 7f are temperature insensitive con- 
stants. See text for details. 



x"(Q + q,^)/XQ = * 



q Twj 



with 



3\2i 



*(x,y) = \x\y/[y z + {x + x A ) 



(21) 



(22) 



It is depicted in Fig. |i"q(b). Although the counter- 
propagating spin wave branches are over-damped by 
electron-hole excitations, they remain two separate 
branches. The intensity is singular at q — and u) = 0, 
and the Lorentzian relaxation energy (the dashed line) 
serves well as the energy scale in this case, as customar- 
ily performed. 



V. SPIN WAVES IN THE AFI PHASE 

While magnetic excitations in the metallic phase are 
dominated by over-damped modes, the AFI phase ex- 
hibits conventional propagating spin wave excitations. 
Word et al. [|l| measured spin waves in this phase around 
(1/2,1/2,0) at 142K in pure V 2 3 . There is a spin ex- 
citation gap of 4.8 meV, which is essentially indepen- 
dent of temperature for T <T N {32). For V1.94Cro.06O3 
(Tat=177K), the gap energy at (1/2,1/2,0) is not differ- 
ent from that seen for the pure sample for temperatures 
below 50K (32). However, it shows substantial softening 
with rising temperature, consistent with the trend to- 
wards a continuous phase transition with increasing Cr 
doping. 
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FIG. 17. Constant energy transfer scans along the (001) 
and (110) directions around the magnetic Bragg point 
(1/2,1/2,0). The two peaks in each frame correspond to 
counter-propagating spin wave modes. The different peak 
shapes are due to resolution effects (focusing and defocusing). 



Here we present spin waves measured at 11 K for 
Vi.944Cr .056O3. Fig. [It] shows constant energy scans 
near a magnetic zone center. The spin wave dispersion 
relation has been measured for energies up to 50 meV 
(see Fig. |l^) which covers about one quarter of the Bril- 
louin zone. For comparison, some of previous data for 
pure V2O3 are also included in Fig. [l^. At low measure- 
ment temperatures and for Tilo < 25 meV, there is no 
difference in the spin wave dispersion for pure and 3% 
Cr doped V2O3. 

Spin waves in the isomorphous corundum antiferro- 
magnets Cr2C>3 and a-Fe^Os have been measured over 
the whole Brillouin zone [107, 



10S ] . These two compounds 



have differen t ma gnetic structures [ 109, 1 10| . Linear spin 
wave theory |mj using the Heisenberg Hamiltonian 



H 



c IT Q z 



(23) 



satisfactorily describe the experimental dispersion rela- 
tion. The dominant exchange constants for Cr2 03 are 
J = —7.5 and —3.3 meV between spin pairs A and B 
respectively [107] (refer to Fig. ||), while for a-Fe203 
spin pairs C and D interact most strongly with exch ange 
constants J = —2.6 and —2.0 meV respecti vely . 1 1 Qg ] 



112 



at 



In V2O3, the corundum symmetry is broken 
the AFI transition. Consequently there are seven differ- 
ent exchange constants between these spin pairs [refer to 
Figfb)]. 

Dispersion relations for spin waves in the AFI phase of 
V 2 03 are derived in Appendix [a| following the procedure 



described by Saenz [113[ . The limited data in Fig. 18 do 
not allow reliable determination of individual exchange 
constants. Along the (hhO) direction the acoustic spin 
wave has a simple dispersion relation 
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FIG. 18. Spin wave dispersion near the antiferromagnetic 
zone center. Solid circles represent the current measurements 
at 11 K for Vi.944Cro.o56 03, the square is from Ref. [ [53]]. 
For comparison, data from [ [53] for pure V2O3 at 142 K are 
represented by open circles. With these limited data, indi- 
vidual exchange constants cannot be derived reliably. Only 
a specific linear combination as described in the text can be 
determined. Insert: the reciprocal lattice plane for the mea- 
surement. The diamonds denote nuclear Bragg points, and 
crosses magnetic Bragg points. 



[H - 4S(Jp H 
-[4S(Jp + J s 



Js + 2J ( )] 2 
f 2J c )] 2 cos 2 (2tt/i). 



(24) 



In combination with the previously measured value of 
the ordered moment (gS = 1.2 ^jj ) our measurements of 
this dispersion relation permits the local anisotropy field, 
H a = 0.13(2) meV, and Jp + J s + 2 J c = 33(2) meV to be 
determined with confidence. [114] To determine the indi- 
vidual exchange constants, experiments which will mea- 
sure spin waves in the AFI across the whole Brillouin 
zone are now under way. 



VI. THE INSULATING SPIN LIQUID 

The solid circles in Fig. [l9| represent constant en- 
ergy scans along (10£) at 205K in the PI phase of 
Vi.944Gro.o5e03. The peak profile is not sensitive to en- 
ergy transfer for 1.6 meV< Huj <18 meV. The peaks 
are very broad and the half-width-at-half-maximum 
(HWHM) corresponds to a correlation length £ c ~ 1.5 A 
in the c direction. Scans in the basal plane (refer to 
Fig. ^) also yielded broad ftw-insensitive peaks, corre- 
sponding to £ a ~ 2.0A in the a direction. These cor- 
relation lengths extend only to the nearest neighbors in 
the respective directions and they are even shorter than 
in the PM phase at comparable temperatures, as can be 
ascertained by comparing Fig. [l^ or Fig. |l5| to Fig. [n]. 
This is a surprising result because the electron-hole pair 
damping mechanism which exists in the metal is inhib- 
ited by the charge excitation gap of the insulator. 

The effect of temperature on the short-range magnetic 
correlations in the PI phase is shown in Fig. EU Data 
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FIG. 19. Constant energy scans along the c direction in 
the PI phase at 205 K (solid circles) and in the AFI phase 
(open circles) for (Vi-xCr^Os (x=0.028). Measurements 
in the AFI phase were conducted at 170 K for hui = 6-12 
meV and at 160 K for other energies. Extra intensity near 
the nuclear Bragg points (102) and (104) for hcu < 2 meV is 
probably due to phonons. Statistical uncertainty in intensity 
is indicated by the vertical bars. The spectrometer configu- 
rations at BT2 were Ei = 12.7 meV with horizontal collima- 
tions 60'-40'-40'-60' for Hoj = 1.63 meV, E f = 13.7 meV with 
60'-40'-40'-60' for Huj = 2-15 meV, and Ei = 35 meV with 
60'-40'-80'-100' for Hlu = 18 meV. 



for two values of energy transfer are shown. Clearly, 
q-dependent magnetic scattering is suppressed at high 
temperatures. Even so, there is no discernible narrow- 
ing of the peak, even as the AFI transition temperature, 
Tjv = 180 K, is approached. Thus, Tn does not ap- 
pear to be the critical temperature for spin fluctuations 
in the PI. As for stoichiometric V2O3 (refer to Fig. fl5| ), 
the characteristic wave vector for magnetic correlations 
of Cr-doped V2O3 moves to another part of reciprocal 
space in the AFI phase leaving a flat background as 
shown by open circles in Fig. [l9| Therefore, the phase 
transition to the AFI is not a conventional antiferromag- 
netic transition, for which the correlation length diverges 
at Tjy and magnetic Bragg peaks emerge from peaks in 
the paramagnetic response function. Instead, the spin 
Hamiltonian is altered by some underlying processes. 

Not only is the location in reciprocal space of the 
paramagnetic neutron scattering in the PI phase un- 
usual, but also is the q and Hlo dependence of scattering 
about this location. For conventional insulating antifer- 
romagnets, when the correlation length is limited to the 
nearest neighbor spin spacing by thermal fluctuations at 
T 3> J/ks, there is very little spectral weight in the mod- 
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FIG. 20. Constant energy scans along the a* direction in 
the PI phase at 205 K for (Vi- x Ct x ) 2 3 (x=0.028). There 
is no appreciable energy dependence of the peak width. The 
spectrometer configurations at BT2 were Ej — 13.7 meV with 
horizontal collimations 60'-40'-40'-60' for hui = 3-15 meV, and 
Ei = 35 meV with 60'-40'-80'-100' for hu = 18 meV. 



ulated part of the dynamic correlation function. In the PI 
phase of Vi.944Cro.o5603, however, constant energy cuts 
through <S(q, ui) show that large spectral weight (> 1A/j,b 
per V, see below) is associated with a ridge along hui 
which forms a very broad peak in q-space. These ex- 
perimental results suggest that the spin-spin interaction 
strength is not dwarfed by the thermal energy and that 
a mechanism other than thermal fluctuations therefore 
must be responsible for limiting magnetic correlations to 
nearest neighbors. 

One possible mechanism is geometric frustration where 
conflicting magnetic interactions prevent long range or- 
der at T r*> J/ks- Solid state chemistry, however, does 
not support this explanation because the dominant ex- 
change interactions in Heisenberg models for c orundum 
antifcrromagnets in general are not frustrated [115, IKS]. 
Disorder introduced by Cr doping is also unlikely to be 
important, since for the same sample at low tempera- 
tures, spin waves in the AFI phase are resolution-limited 
(refer to section |v| and Fig. 3 in Ref. |37| ) . A very promis- 
ing explanation, based on the electronic state of V2O3 
i n th e single ion limit, was recently proposed by Rice 
[117,B7J. Each magnetically active electron on a vana- 



dium ion can exist in one of the two doubly degenerate 3d 
orbitals. The orbital degeneracy may be the sou rce o f 
large orbital contribution to the nuclear relaxation [118|. 
Depending on the relative occupation of these two de- 
generate orbitals on neighboring sites, the exchange in- 
teraction between spins on these sites can have different 
signs [HE]. In the PI phase, the orbital occupation is 
fluctuating and disordered, resulting in a spin Hamilto- 
nian with fluctuating exchange interactions of different 
signs which induce local antiferromagnetic correlations, 
while dephase longer rang e correlations. We will return 
to this idea in section VTJ. 
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FIG. 21. Temperature dependence of constant energy scans 
in the PI phase of Vi.944Cro.o5603. The left panel is for 
hui — 3 meV, and the right panel for hui = 15 meV. The spec- 
trometer configurations at BT2 were Ef = 13.7 meV with 
horizontal collimations 60'-40'-40'-60'. The solid line is a fit 
to the product of the form factor |/(q)| 2 and a gaussian func- 
tion. 



A. Analysis of spatial correlations 

A consequence of the fcu-insensitive peak width in the 
PI phase (Fig. 



19 



and |2fj) is that constant energy scans di- 
rectly probe the q-dependence of the instantaneous spin 
correlation function. Because of the short correlation 
lengths, the spatial dependence of the magnetic corre- 
lations can be modeled by a finite spin cluster which in- 
cludes the four types of near neighbor spin pairs: A, £?, C 
and D [refer to Fig. |3|(a)]. In other words, the summation 
in the structure factor 



5(q) 



{g/j, B ) 2 n 



(S ) 2 +5>os(q-R')(So-S R , 



R' 



(25) 



needs only include one A-type spin pair, three -B-type 
spin pairs, three C-type spin pairs and six D-type spin 
pairs for any one spin, So, of the n = 12 spins in a 
hexagonal unit cell. 

In Fig. |2^, an extended scan along the c axis, covering 
more than three Brillouin zones, and scans in the basal 
plane intersecting the peaks near the (101) and (105) 
are shown for Vi. 944(^0. 056O3 at 205 K in the PI phase. 
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FIG. 22. Constant-?^ = 3 meV scans (a) along (10/), (b) 

along (ft, 0,1), and (c) along (/i,0,5) at 205 K in the PI phase 
(solid circles) and at 160 K in the AFI phase (open circles). 
Squares in (c) represent data at 400 K. (1,0,8), (1,0,2), (1,0,4) 
and (1,0,10) are nuclear Bragg points. The experimental con- 
figuration is indicated in frame (a). 



Background measured at 160K in the AFI phase is also 
shown. A good approximation to the magnetic intensity 
at 205 K is obtained by subtracting the 160 K background 
and such difference data are shown in Fig. 23. The dotted 
line indicates the magnetic form factor for the V 3+ ion 
pEfl , |/(q)| 2 - Extra intensity at large |q| is attributed to 
incomplete subtraction of phonon neutron scattering. A 
simultaneous fit of |/(q)| 2 S(q) to the magnetic intensity 
from these 3 meV scans is shown as solid lines in Fig. 23. 
As can be seen in the figure this spin cluster model is 
consistent with the data over the large part of q-space 
where the magnetic form factor is appreciable. 

The near-neighbor spin correlations derived from the 
best fit (x 2 = 1.9) are: 

(So • S [0 ,o,i/6+5]>"V< S } 2 = 0.08(9), (26a) 

(So ' S [1/3 ,W B /<S) 2 = -0-10(4), (26b) 

(So • S [2/ 3 4 /3, f 5-i/ 6 ]) C /(S> 2 - 0.04(4), (26c) 

(So • S [2/ 3, 1/3 ,i/ 6 ]) D /<S) 2 = -0.06(2), (26d) 

where S — 0.026 and superscripts A-D refer to spin pairs 
in Fig. ||a). Now, it is likely that warming yields an 
additional non-magnetic background, most probably due 
to phonon and multiphonon scattering responsible for the 
upturn at high |2| in Fig. ^2](a). We account for this 
rising background by incorporating an additive constant 
in our fits, as we have done previously [ j35T]37jl . The best 
least squared fit (x 2 = 1.2) is obtained with magnetic 
parameters |35|j37| 

(So • S [0A1/6+5] ) A /(S) 2 = 0.6(3), (27a) 
(So • S [lA2/3j] ) B /(S) 2 - -0.19(8), (27b) 
(So • S [2AlA5 -i/ 6] ) c /<S) 2 = 0.18(8), (27c) 
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FIG. 23. Magnetic signal at 205 K in the PI phase, af- 
ter subtracting the smoothed 160 K intensity (solid line in 
Fig. ^) as a measure of non-magnetic backround. The back- 
ground for (c) was determined from the 400 K and 160 K data. 
The dotted line represents the magnetic form factor |/(q)| 2 . 
Extra intensity at large |q| is probably due to phonon scat- 
tering. The solid line represents |/(q)| 2 <5(q) calculated with 
the parameters of Eq. (26) while the dashed line is a constant 
plus this quantity calculated with the parameters of Eq. (p7[). 



(S -S 



[2/3,1/3,1/6]/ 



-0.09(3), 



(27d) 



and a constant background of 25 counts in 10 min. 

The physical picture which emerges from both fits is 
that in the PI phase, the nearest neighbor spin pairs in 
the basal plane, which contains a puckered honeycomb 
lattice for V, are more likely to be anti-parallel to each 
other, while the nearest neighbor spin pairs in the c di- 
rection are more likely to be parallel to each other. While 
the signs of correlations for different pairs are clear, the 
specific values of the spin pair correlations have large sta- 
tistical error bars and possibly systematic errors associ- 
ated with uncertainty in the background determination. 
If this dynamic short-range spin correlations were allowed 
to develop into a long-range order, the spin correlations 
would be 

(So ■ S[ ,o,i/6+5]) A / (S) 2 = 1, 
(So • S [1/3i2/3 j ] ) B /(S) 2 = -1, 



(So • S[2/3,l/3,5_l/6]) / (S) 2 
(S ' S[ 2 /3,l/3 > l/6])' D / (S) 2 



-1. 



This corresponds to the type I collincar antifcrromagnetic 



structure for corundum [llf 

In Fig. ^J, the magnetic intensities of Fig. |19| from 
hui = 1.63 meV to 18 meV have been integrated, to yield 
the quasi-instantaneous correlation function 

(.9MB ) 2 ■ 



h / c/w5(q, uj) 



3N 



RR' 



q ' (R '- R) (S R (ft)-SR-(Q)) 
(28) 
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FIG. 24. Energy integrated intensity, up to 18 meV, as 
a function of q = (W£) at 205K in the PI phase for 
Vi.944Cro.o5603. The solid line is 125(q) based on Eq 
the parameters of Eq. (26), and s 2 S 2 /3 = 0.80 2 . 
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FIG. 25. (a) The dynamic susceptibility at Q specified in 
the figure, and (b) the local dynamic susceptibility as a func- 
tion of energy for Vi.944Cro.o5603 at 205K in the PI phase. 
The corresponding quantities for stoichiometric V2O3 at 200K 
in the PM phase are shown in (c) and (d). 



where \6t\ < n/18meV= 3.7 x lCr 14 sec. The solid line in 



the figure is <S(q), calculated from Eq. (|25|), (26) using 
(5Mb) 2 (S> 2 /3 = (0.80^b) 2 . Thus, at 205K, the fluctuat- 
ing moment involved in the short range correlations, up 
to hu> = 18 meV, is 1.4(2)/iB per V. This is as large as 
the ordered moment of 1.2/xs per V in the AFI phase. 



B. Energy dependence 

Fig. ^||(a) shows the energy dependence of x"(Q> w ) a t 
Q = (1,0,0.76) for Vi. 9 44Cr . 056 O 3 at 205 K. The data 
were obtained from the peak intensities of the constant- 
energy scans in Fig. [nj. There is a remarkable similarity 
between these data and the peak intensities for metallic 
pure V2O3 shown in Fig. |25|(c) despite the difference in 
5(q) for these two samples (see, e.g., Fig. [l2| and |l9|) . 

The spectrum can be modeled by a Lorentzian 



x"(q,w) = 



x(q) 



Tl/jJ 



r q i + (^/r q ) 2 ' 



(29) 



where x(q) is a generalized q-dependent static suscep- 
tibility and r q is the relaxation energy. The data in 
Fig. |25|(a) only allow us to place a lower limit of 18 meV 
on r Q at T = 205 K. The initial slope 

II(Q) = lim X "(Q,w)/^ = x(Q)/r Q (30) 

is however reliably determined to be n(Q) = 
0.075(4)/i|/meV 2 per unit cell. 

The energy dependence of the local magnetic suscep- 
tibility 



X"(w) = J d 3 q X "(q,c)/ J d 3 q, 



(31) 



for the PI and PM are shown in Fig. |25|(b) and (d). For 
the insulating sample these data were derived from the 
Gaussian fits shown in Figs. |l9| and ^ while we used the 
fits shown in Fig. |l5| for the metallic sample. Again we 
observe that for T 200 K the magnetic spectral weight 
for hu) < 25 meV is substantially larger in the PI phase 
than in the PM phase. 



C. Comparison to bulk susceptibility data 

As a check on the reliability of our normalization pro- 
cedures, we compare the normalized dynamic suscepti- 
bilities at T — 205 K discussed above to the conventional 
bulk susceptibility, x = x'(0>0) = x(l = 0) + x'(0, 00) 
for Vi.944Cr .056O 3 . 

From Eq. ( |29| ) and the Kramers-Kronig relation, as 
long as we are in the classic regime, 



r q < k B T, 



S(q) = fi J du X "(<l,u) + 
Therefore, at a given temperature, 



fc B Tx(q). (32) 



x(o) 



g(g) 
5(Q) 



x(Q) = 



S(0) 
S(Q) 



r Q n(Q). 



(33) 



Using parameters in (26) or ( |27| ) for our finite cluster 
model, we find 5(Q)/5(0) w 3.0 or 2.2. Thus, assuming 
that Tq < kgT, we can estimate the bulk susceptibility 



x(0) from the parameters quoted in section VI B. The 
result is x > 0.5-0. 7/i 2 3 /meV per unit cell = 1.1-1.6 x 
10~ 3 emu/mole V for Vi^Cro.oseOa at T = 205 K. 
This number should be compared to the measured bulk 
susceptibility which is 3.6 x 10 -3 emu/mole V at 205K. 

M 



D. Temperature dependence 

The temperature dependence of x"(Q> i*>)/7ia; is shown 
in Fig. Effl for Vi.g44Cro.o5603- The 3 meV and 15 meV 
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FIG. 26. Temperature variation of x"(Qi at 
Q = (1,0,0.76) for Vi.944Cro.o5603 in the PI phase, with 
Tiw used in the measurement indicated in the figure. The cor- 
responding quantity at the low energy limit, Xq/jak 2 , for 
V1.973O3 in PM (diamonds) is included for reference. The 
solid line is proportional to T -2 , the dotted line is propor- 
tional to (T 2 — (180 if) 2 )" 1 , and the dashed line is propor- 
tional to (T 2 — (9 AT) 2 ) -1 . 



data are from Fig. |2l] and the 9 meV data from Fig. 4(c) 
of Ref. j37|. The values of x"(Q,w)/&w from different 
hui = 3-15 meV are close to each other down to Tjv = 
180 K. This indicates that x"(Q,w) cx hui in this energy 
range and that Tq therefore exceeds 15 meV even close 
to the AFI transition (refer to Eq. (|29|)). This also means 
that the ratios plotted are a measure of II(Q) [Eq. (fjOj)]. 

In a second-order antiferromagnetic transition, x(q) cx 
{T+T^lL 1 an d Tq oc T—Tpf in the Gaussian approxima- 
tion pfjj , therefore, n(Q) cx (T 2 — T^ ] 2 )~ 1 scales as T~ 2 



for T ^> T^r. As expected, experimental data for metallic 
V1.973O3 follow this relation for T > T N = 9 K (refer 
to diamonds in Fig. ||^). Interestingly, the data for the 
temperature dependence of II(Q) in Vi.g44Cr .056O3, are 
proportional to T~ 2 (refer to the solid line) all the way 
down to Tpf — 180 K. Therefore, the spin system in the 
PI phase at temperatures even slightly above the PI- AFI 
phase transition appears to be far away from the critical- 
ity associated with the real T/y. On the other hand, 
is remarkably close to Neel temperature of the metallic 
AFM found for hole-doped samples. In other words, the 
PI fluctuations are those which one would associate with 
nearly the same quantum critical point as found for the 
PM and AFM phases. 



VII. DISCUSSION AND CONCLUSIONS 

In the previous sections, the magnetic properties of 
V2O3 in each of its phases were described separately. In- 
terest in this material, however, is driven by attempts to 
understand the Mott-Hubbard transition. Here, we will 



discuss the implications of our magnetic neutron scatter- 
ing experiments for the interpretation of the transitions 
between the PI, AFI, PM, and AFM phases. 

The discovery of antiferromagnetic order in metallic 
V2O3 through the 57 Fe Mossbauer effect prompted 
much theoretical work on antiferromagnetism in strongly 
correlated Fermi liquids. Phase diagrams including 
a metallic antiferromagnetic phase for the Mott sys- 
tem were subsequently published by several authors 
p8| ^,^2). However, for many of these studies, a Heisen- 
berg antiferromagnetic state like that in the insulator 
was put into the metallic state by hand and local mo- 
ment spin dynamics in this phase was implicit. Even 
though the effective mass of electrons in metallic V2O3 
is strongly enhanced, i.e., the Fermi liquid is nearly lo- 
calized near the metal-insulator transition, we showed 
in the previous sections that the antiferromagnetism in 
metallic V2O3 is completely different from the antiferro- 
magnetism in the AFI. While Heisenberg models of ex- 
change coupled unpaired spins located at V ions can ac- 
count for magnetic phenomena in the AFI phase, antifer- 
romagnetism in the metal is controlled fundamentally by 
the Fermi sea, j3^]3(|[7Cj . Specifically the small-moment 
static SDW order results from a nesting or near-nesting 
Fermi surface, and the broad bandwidth spin excitation 
spectrum is reminiscent of an electron-hole pair contin- 
uum. In short, our experiments show that the metal- 
insulator transition separates two qualitatively different 
antiferromagnetic states, with localized moments in the 
AFI and itinerant moments in the AFM phase. 

With the large number of spin fluctuation modes in 
the Stoner continuum, spin fluctuations in an itinerant 
antifcrromagnet have a stronger renormalizing effect on 
equilibrium magnetic properties than spin waves have in 
an insulating antiferromagnet. This makes magnetism 
at finite temperatures for itinerant magnets much more 
complex than for localized moment magnets. Near the 
quantum critical point in the small moment limit, scaling 
relations can be derived based on general physical argu- 
ments S-IbtJ. The scaling functions and relations be- 
tween non-universal constants can be further calculated 
by the SCR theory @||Jl|JlO|. We found remarkable 
quantitative consistence between experimental results for 
metallic V2-j,C>3 and the theory. 

The source of an itinerant antiferromagnetic insta- 
bility, which is characterized by singularity of general- 
ized magnetic susceptibility at finite q, could be a di- 
vergent Lindhard function due to nesting Fermi surface 
]67| , p9]| , a nearly nesting Fermi surface in combination 
with the Coulomb interactions {39, 121 1, or solely the 
strong Coulomb interactions. The last possibility is still 
under theoretic investigations. [122 123 The nesting sce- 
nario has been extremely successful in describing the clas- 
sic itinerant antiferromagnet, Cr. [n| We find it also 
pro vides a comprehensive framework for metallic V2O3 
P3|)p6[ . In contrast to Cr, however, only a small area of 
the Fermi surface is involved in the SDW Q , and strong 
correlations in metallic V2O3 require only near nesting 
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of the Fermi surface for the long-range antiferromagnetic 
state to occur. 

Three kinds of metal-insulator transitions occur in the 
V 2 3 system, the PM-PI transition, the PM-AFI transi- 
tion, and the T^O AFM-AFI transition. Whether they 
should be called Mott transitions has been controversial 
and, to some extent, is confused by the evolving mean- 
ing of the Mott tran sition. While previously noticed 
by others, Mott 124] was responsible for bringing at- 
tention to the possibility of localization of electrons in 
a narrow band through electronic correlations, which is 
neglected in conventional band theories of solids. Mott 
argued for a first-order metal-insulator transition when 
the electron density is reduced below some critical value. 
Not much detail about the electronic processes taking 
place at the transition were known at th at ti me and since 
the "exceptional" insulators like NiO |l25| , which mo- 
tivated Mott's research and are today called Mott and 
charge-transfer insulators, were known experimentally to 
be antiferromagnets, it was unclear whether there was 
any new physics apart from magnetism. Specifically, 
if these insulators are always antiferromagnetic, a dou- 
bling of the unit cell will open a gap at the Fermi level 
as shown by Slater |27|. However, Hubbard ||] demon- 
strated in 1964 that a metal-insulator transition can be 
produced by strong Coulomb correlations without anti- 
ferromagnetic ordering. Brinkman and Rice ^] in 1970 
outlined a distinct type of Fermi liquid on the metal- 
lic side of the Mott-Hubbard metal-insulator transition. 
Strong electron correlations in the PM phase of V2O3 
were revealed through its Brinkman-Rice-like behavior 
[126 1 and the PM-PI transition was discovered as a pos- 
sible experimental realization of the M ott-H ubbard tran- 
sition by McWhan, Rice and Remeika [ 127 1 . This proba- 
bly stands as the first explicitly identified Mott-Hubbard 
transition in a real material. 

The difference between the paramagnetic metal and 
the paramagnetic insulator in (Vi-^Cr^Os disappears 
at the critical point near 400 K ( or th e critical line in 
the composition-P-T phase space [128]). A continuous 
crossover between paramagnetic spin fluctuations in the 
metal and the insulator thus is expected on a path above 
the critical point. By inspecting the evolution of mag- 
netic correlations with temperature and doping in the 



metal (refer to Figs. |l 
tor (refer to Figs. [l9] and 



and |l5| ) and in the insula- 
21]), such a magnetic crossover 
from metal to insulator can be conceived as follows: The 
pair of incommensurate peaks in the metal broaden in 
q upon approaching the critical point until they finally 
merge into the single broad peak which we observe in the 
PI. The blurring of the incommensurate peaks which is a 
measure of Fermi surface dimension on the nesting parts 
parallels the losing of metallic coherence among electrons. 

Although Hubbard as well as Brinkman and Rice have 
shown that a Mott transition does not require an an- 
tiferromagnetic transition, antiferromagnetic order does 
promote the insulating state. This is apparent because 
of the enlarged insulating phase space associated with 



the AFI in various studies based on the one-band Hub- 
bard model |2^ , |29| , |l2]| . The resemblance of the theoreti- 
cal phase diagrams to the experimental phase diagram of 
V2O3 gives merit to Slater's idea of an insulator induced 
by antiferromagnetism. 

However, an important difference exists between these 
one-band phase diagrams and the V2O3 phase diagram. 
The magnetic transition in the theories is predicted to be 
of second order while the experimental AFI transition in 
V2O3 is first order [refer, e.g., to Fig. ||(a)]. This differ- 
ence is not trivial such as being merely a magnetostric- 
tion effect. Beneath it is a fundamental difference in 
temperature depe nden ce of spatial spin correlations |37| . 
It is well known [ |129| that a one-band Hubbard model 
reduces to a Heisenberg model with exchange J ~ t 2 jU 
in the insulating limit U/t 3> 1. For temperatures lower 
than the insulating gap ~ U/ks, the low energy physics 
is contained in the Heisenberg model with a second or- 
der antiferromagnetic transition at T/v ~ J/fc^. The 
characteristic magnetic wave vector remains the same for 
T > T N and T < T N . In V 2 3 , dynamic spin fluc- 
tuations in the corundum PM and PI phases peak at 
a magnetic wave vector along the c* axis. These mag- 
netic correlations abruptly vanish and are replaced by 
(l/2,l/2,0)-type long range order and the associated spin 
waves below the first-order AFI transition (see Ref. [ [B7J ] 
for details; see also Figs, [l^ and [H] in this paper). This 
switching between spin correlations with different mag- 
netic wave vectors indicates that the antiferromagnetic 
transition is not a common order-disorder phase transi- 
tion, as is found, for example, in the parent compounds 
of high-Tc superconductors, but that instead the spin 
Hamiltonian is somehow modified at the first-order tran- 
sition to the AFI phase. These results find no compre- 
hensive explanation in one-band Hubbard models. 

We mentioned in section VI that the electronic state 
of V2O3 in the single-ion limit can be approximated by 
two degenerate d-orbitals at each V ion site which are 
filled by a single electron. New physics in Hubbard mod- 
els with degenerate bands was di scus s ed by se veral au- 
thors more than twenty years ago [13C-132,119[. It gains 
renewed inte rest in light of recent experimental progress 
[ 117,133 , 134 1 . Depending on which one of the degenerate 
orbitals is occupied at a given site and at its neighboring 
sites, the exchange interaction between the nearest neigh- 
bor spins is different and can even change sign. Thus, 
besides the spin degrees of freedom, a new orbital de- 
gree of freedom enters the low energy physics. The spin 
and orbital degrees of freedom are strongly coupled, and 
the effective spin Hamiltonian depend s sensiti vely on or- 
ders in the orbital degrees of freedom. [ 117, 1 19 1 Switching 
between spin correlations with different magnetic wave 
vectors was also recently observed in manganite at an 
orbital/charge ordering transition. 135 

The anomalously short-range spin fluctuations in the 
PI and the abrupt switching of magnetic correlations at 
the AFI tran sitio n can be understood, as first pointed 
out by Rice 117], if the primary order parameter for 



19 



the AFI transition involves orbital degrees of freedom. 
At low temperatures, order in orbital occupations breaks 
the three-fold corundum symmetry. Pairs with differ- 
ent orbitals at neighboring sites produce different bond 
strengths and different exchange interactions, resulting in 
the structural distortion and the antiferromagnetic struc- 
ture in Fig. U (b) which has two antiferromagneti c bo nds 
and one ferromagnetic bond in the basal plane. [119| In 
the PI phase the orbital occupation fluctuates between 
the two states of the doublet, thus restoring the corun- 
dum symmetry and producing a different spin Hamilto- 
nian with random and fluctuating exchange interactions. 
With this kind of spin Hamiltonian, magnetic correla- 
tions necessarily are limited to the nearest neighbors. 
The fact that the magnetic transition is strongly first or- 
der indicates that the intrinsic Neel temperature, T^ FI , 
of the spin Hamiltonian associated with the AFI phase is 
larger than the orbital ordering temperature To- Thus, 
when orbitals order and the spin Hamiltonian changes 
from the one associated with the PI phase, which has a 
T^ 1 < Tq, to the one associated with the AFI phase, 
magnetic moment reaches from zero directly to a value 
appropriate for T = T < T$ FI ■ 

In conclusion, we have shown that the metal-insulator 
transition fundamentally changes the nature of antiferro- 
magnetism in the V2O3 system. Metallic V2O3 is a pro- 
totypical strongly correlated small-moment itinerant an- 
tiferromagnet, involving a nesting instability on a small 
part of Fermi surface. We also have presented strong indi- 
cations that orbital degeneracy plays a crucial role in the 
magnetism and phase transitions of this material. Specif- 
ically it appears that these orbital degrees of freedom give 
rise to a novel fluctuating paramagnetic insulator, and a 
first order transition to an antiferromagnetic phase with 
a spin structure which is unrelated to the paramagnetic 
short range order. At the same time, the fluctuations 
in the paramagnetic insulator and paramagnetic metal 
are so similar that they suggest that the electron-hole 
pair excitations associated with Fermi surface nesting in 
the metal are transformed into orbital fluctuations in the 
insulator. 
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APPENDIX A: SPIN WAVES IN AFI 

Spin wave excitations from a collinear spin arrange- 
ment were studied by Saenz [US]. For a material with 
n magnetic ions in a primitive magnetic unit cell, the 
positions of these ions are described by 



R 



I ft 



R, 



(Al) 



where Ri is the location of the ith unit cell (i=I,2,...,N), 
and 17, is the position of the /nth ion inside the cell 
(fx = 1,2, ...,n). The vector spin operator for the ion 
denoted by {i/i} is Si M and the spin quantum number is 
Sfj,. Choosing z to be parallel to the collinear direction, 
the spin orientations are described by e M = ±1 for up or 
down spins in a unit cell. A diagonal n x n matrix, E, 
with elements e M 6 M „ is introduced for later use. 
The spin Hamiltonian is given by 



T~t — Jifiji/Sin • y] e^H^Sf^, (A2) 



where models the anisotropy field and the collinear 
component of the external field on the /ith ion, and the 
exchange interaction Ji^.j V depends only on fx, v and the 
ionic separation Ri M — Hj V . As usual, Ji^.jy — Jjv,in and 

The lattice transformation of the exchange constants 
is a n x n Hermitian matrix function and its elements are 

■V ( K ) = X! J o^fin expfiK • (Rjv - R 0m )], (A3) 



and k is defined within the first Brillouin zone. Introduce 
the matrix L(re): 

L^( K )42(e M ^J^(0)e„5, J + \h^5^ 



(A4) 



which is Hermitian V(k) — L(k), and L(— k) = L*(/c). 
The n eigenvalues, A M (k), of EL are related to the energies 
of the n spin wave modes by 



hw^K) = |A jLI (k) 



(A5) 



For a stable spin structure, the matrix L(/c) is non- 
negative definite for every k, and A^(k) does not cross 
the zero axis for any branch of the dispersion relation, 
i.e. it is either non-negative or non-positive throughout 
the Brillouin zone. The number of non-negative (non- 
positive) branches of A M (n) equals the number of positive 
(negative) elements of e^. 
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The spin Hamiltonian now can be rewritten in terms 
of spin wave modes 



H = -N 



(A6) 



in the semi-classic limit [136], namely, — £n(S?) <C 

The primitive magnetic unit celljplj in the monoclinic 
AFI phase of V2O3 is shown in Fig. ||(bj together with the 
conventional hexagonal cell of R3c symmetry. Neglecting 
the small lattice distortion in the phase transition to the 
AFI, the monoclinic and hexagonal unit cells are related 
by the transformation 



(A7) 
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-1/2 1/2 




(A8) 
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-2 / 


c™ 





Denoting the above two matrices as A and B respectively, 
the corresponding transformation matrices for reciprocal 
lattice primitive vectors are B T and A T . 

There are four magnetic V ions per monoclinic unit 
cell [refer to Fig. ||(b)] and they are assumed to have 
equivalent spins, 5 M = S. With the designation of ions, 
1-4, in Fig. |](b), (e^) = (1, 1, —1, —1) for the experimen- 
tally found magnetic structure in the AFI. The local field 
(Hp) = (H+,H + ,H-,H-) with H ± = H ± H e , where 
Hq is the anisotropic field in the material and H e is an ex- 
ternal field applied along the up spin direction. Limiting 
the range of exchange interactions to the fourth nearest 
neighbors, there are seven distinct exchange constants: 
J a , Jp, J 7 , Js, J c , J( and J v (refer to Fig. ||(b)). Their 
contributions to various spin pairs are tabulated in Ta- 
ble 0. Now the matrix J (k) can be readily calculated by 
collecting contributions to spin pairs ((iv) according to 
([A3]). Notice that spin pairs (/z/i) are identically con- 
nected by J,,, therefore Jn=J22=J33=J44- Spin pairs 
(23), (24) and (34) are identically connected as spin pairs 
(14), (31) and (21) respectively. This leads to J23=Ji4, 
J24=J3i and J34=J2i- Recall also that J(k) is Hermitian, 
the matrix has only four independent elements: 



J(«) 



They are given by 



A B C D 
B* A DC* 
C* D* A B* 
D* C B A 



(A9) 



TABLE II. Spin pairs connected by the seven exchange 


constants. 
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1 z x 




6 4 —X 
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toy — c m /z 




1 3 y + c m /2 




2 4 — y — c m /2 




2 4 — y + c, n /2 




1 2 X 3 777, 


Js 


1 3 y — a m — c m /2 




1 3 y — a m + c m /2 




2 4 a m — y — c m /2 




2 4 a m — y + c m /2 
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1 z x + b m 


7 

■k 


1 4 b m /2 — c m /2 




1 4 b m /2 + c m /2 




1 4 —Om/2 — C m /Z 




1 4 — b m /2+Cm/2 




2 3 b m /2 — Cm/2 




2 3 b m /2 + c m /2 




2 3 —bm/2 — c m /2 




2 3 —bm/2 + c m /2 




1 1 a m 




1 l c -a m 


a x = R02 — Roi 


= (l/3 + 2<S)(a m -b m ), where 5 = 0.026. For 


M=(21), R,„ 


— Ro M = —x. Exchanging fi and v switches 


the sign for Rj„ 


— R()m- 


b y = (1/3 + 25) 


a m + (1/6 - 2(5)b m . 


C J,, has identical contributions for (jt«/)=(22), (33) and (44) 


spin pairs. 




A = Jh(k) = 


2J v cos(2irh), (AlOa) 


B = Ji 2 (k) = 


J Q exp[i27T<^>(/i — k)] + 




J 1 exp[z27r(— iph — (pk)] + 




J e exp[i2n((f)h + ipk)}, (AlOb) 


C = Jl3(#6) = 


2 Jp cos(tt/) exp[z27r(0/i + 9k)] + 




2J S cos(ttZ) exp[i2?r(-'0/i + 0fc)], (AlOc) 


D = J 14 (k) = 


4J C cos(7rZ) cos(7rfc), (AlOd) 


where k = (hkl) is indexed in the monoclinic reciprocal 


lattice, <f> = 26 


+ 1/3, %jj = 2/3 -25, 9 = 1/6 - 2(5 and 


6 = 0.026. Notice that only B and C are complex. A 


and D are real 




The diagonal terms of L in the first line of (|A4j) are 


6^[H^ + 2S(J a + J 7 + J e + 2J r) - 2J P - 2J S - 4J f )]. 


When there is 


no external field, Hn = Hq no longer de- 


pends on [i. In this case, 



L(k) = 25 
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(All) 



where 
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F = J a + J 7 + J e + 2 — 2 Jg — 2 



The eigen value equation for EL is 

F — X' —B C D 

-B* F — X' D C* 

C* D F + X' -B* 

D C -B F + X' 



4J C + H /2S - A. 

(A12) 



(A13) 



where A' = A/25. It can be verified that the linear and 
cubic terms in A' vanish and therefore the eigen value 
equation is quadratic in A' 2 : 

A' 4 -2b( K )X' 2 + |L(k)|=0, (A14) 

where both the determinant |L(k)| and 

6(#c) = F 2 + \B\ 2 - \C\ 2 - D 2 (A15) 

are real functions of k. The consequence of this is that in 
the absence of an external magnetic field, there are two 
branches of doubly degenerate spin wave modes with a 
dispersion relation 

Huj(k) = 2S[b(n) ± c(k)] 1 / 2 , (A16) 

where 

c 2 (n)^b 2 (n)-\L(K)\ 

= AF 2 \B\ 2 + A\C\ 2 D 2 - 2\B\ 2 \C\ 2 + 

B 2 C* 2 + B* 2 C 2 + AFD(BC* + B*C). (A17) 
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